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\0 ' ABSTRACT 



. Context. X-ray studies of active galactic nuclei (AGN) with powerful nuclear winds are important for constraining the physics of the 

' inner accretion/ejection flow around supermassive black holes (SMBHs) and for understanding the impact of such winds on the AGN 

environment. 

Aims. Our main scientific goal is to constrain the properties of the circum-nuclear matter close to the SMBH in the mini-broad 
. absorption line quasar (mini-BAL QSO) PG 1 126-041 using a multi-epoch observational campaign with XMM-A'ewton. 

Methods. We performed temporally resolved X-ray spectroscopy and simultaneous UV and X-ray photometry on the most complete 
set of observations and on the deepest X-ray exposure of a mini-BAL QSO ever. 

Results. We found complex X-ray spectral variability on time scales of both months and hours, which is best reproduced by means 
of variable massive ionized absorbers along the line of sight. As a consequence, the observed optical-to-X-ray spectral index is found 
^ ' to be variable with time. In the highest signal-to-noise observation we detected highly ionized X-ray absorbing material outflowing 

' much faster {vx ~ 16 500 km s"') than the UV absorbing one (f,,,, ~ 5 000 km s"'). This highly ionized absorber is found to be 

variable on very short (a few kiloseconds) time scales. 
■ Conclusions. Our findings are qualitatively consistent with line-driven accretion disk winds scenarios. Our observations have opened 

^ ' the time-resolved X-ray spectral analysis field for mini-BAL QSOs. Only with future deep studies will we be able to map the dynamics 

\Q ' of the inner flow and understand the physics of AGN winds and their impact on the environment. 
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2 1. Introduction dHaardt & Mara schil [1991 iHaardt et al.l 11994 . What is clear 

I ^, r , ■ ■ / , , ^ n ■ both from variability and microlensing studies is that the X-ray 

The structure of the mner regions (sub-parsec scale) of active ^^.^^.^^ .^^ ^^^^ ^^^^^^^ ^^^^ on e, and it spans 



X 



galactic nuclei (AGN), as probed by Uy and x-ray observations pniy a few up to tens of gravitationa l radii (e.g.. iMorgan et all 
seems to be very complex, and certamly it is still not understood io^; Chartas et al. 2009b; Dai elMMM- An extensive X-ray 



well. The optical/UV continuum emission is most probably due ^^^^^^^ inner regions of the nearby Seyfert 1 .8 NGC 
^ to the thermal emiss ion from an optical y flnc k, geometrically 1365 has made it possible to constrain the size of the X-ray emit- 
. - . thin accretion disk ^Shakura & Sunvaeyl 092^. The disk sur- ^. .^^ ^^^.^ ^ ^ ^^^^jj independent tech- 
rounds the central supermassive black hole (SMBH), spanning ^. , 2009b). The UV and X-ray continuum pho- 
radii from a few up to several hundreds of gravitationa radii ^^^^ ^^^^^^^^ reprocessed (absorbed, re-emitted, scat- 
(r = GMbh/c^), I.e. possibly from the innermost stable cir- ^^^^^^ reflected) by gas and dust in the inner regions of AGN, 
culai- orbit ai-ound the SMBH up to the disk self-fragmentation ^^^^ ^j^^^^ ^^^^^ ^^^^^^^ -^^^ continuum 
radius (e.g., Prm gle„ 1981 ) The origin of the X-ray continuum ^^.^^.^^ ^^^^^ ^^.^^^^ ^^^^ ^ ^ comprehensive analysis 
emission IS less clearly understood and is thought to be the re- 2005b a, 2007, 2009a, again on NGC 1365). 
suit of the Comptonization of accretion disk UV seed photons Conversely, observing the reprocessing features in the UV and 
into a cloud of very hot electrons, the so-called X-ray corona ^-ray bands can put strong contraints on the geometry, the phy s- 
Send offprint requests to: M. Giustini, e-mail: i^al characteristics, and the dynamics of the gas hosted in the 
giustini@iasfbo . inaf . it inner regions of AGN. In the past few decades, spectroscopy in 
* Based on observations obtained with XMM-Newton, an ESA sci- die UV and X-ray band have revealed the presence of substan- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA. 
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tial column densities of ionized gas flowing out from the inner 
regions of AGN. 

In the UV band we observe blueshifted spectral absorption 
lines due to resonant transitions of ionized metals such as Mg II, 
Al III, Si IV, C IV, N V, O VI. Depending on the width of 
the absorption troughs, quasars hosting such features are clas- 
sified as broad absorpti on line quasars (BAL QS Os, FWHM 
> 2000 km s ', e g. iTurnshek et al.1 [1980; We vmann et al.l 
ll99lUReichardetal.ll2003h . mini-broad absorption line quasars 
(mini-BAL QSOs, 500 km s"' < FWHM < 2000 km s"'), 
and narrow absorptio n line quasars (NAL QSOs, FWHM 
< 500 km s '; see 'Gangul y et all 1200 it iHamann & Sabral 
12004) . BALs features are observed in ~ 15% of optically se- 
lected QSOs (e.g.. iHewett & Foltzl 1200 31; iReichard etal] 120031; 
iKnigge et al.' '2008'; 'A llen et al.1 1201 ih . Mini-BALs and NALs 
are together observe d in ~ 12 - 3 0% of optically selected 
QSOs (Gang uly & BrothertonI |2008|) . Despite the very differ- 
ent absorption trough widths observed, BAL, mini-BAL, and 
NAL QSOs share the same velocity ranges, reaching UV ter- 
minal velocities vl'J, from a few 10^ km s"' up to several 
10'^ kr n s~' (e.g., Rodriguez Hidalgo et al. 200_7; Misawa et al.. 
l2007at iGibson et a l. 2009b). Furthermore, is known to cor- 
relate with the continuum lum inosity Luy ( iLaor & Brandtl200l 
iGangulv & Brothertonll2008l) . meaning that the highest terminal 
velocities are observed in the objects with the highest UV con- 
tinuum luminosity, suggesting an important role of the AGN UV 
radiation pressure in the acceleration of such winds. 

In the X-ray band absorption due to ionized species, such as 
N VI- VII, O VII- VIII, Mg XI-XII, Al XII-XIII, Si XIII-XVI, as 
well as L-shell transitions of Fe XVII-XXIV, are observed to be 
blueshifted by a few hundred to a few thousand km s~^ in ~ 50% 
of type 1 AGN (the "warm absorber' ', e.g. iRevnoldsl [19971 
iPiconcelli et all 120051: iMcKeman et al.ll2007h . Moreover, in re- 
cent years, thanks to the high collecting area of X-ray satellites 
such as XMM-Afewfon, Chandra, and Suzaku, blueshifted ab- 
sorption lines due to highly ionized gas (i.e. Fe XXV, Fe XXVI) 
outflowing at much higher velocity {vout ~0. 05-0. 2c) have been 
observed in a number of AGN, at both low and high redshift (see 
e.g. Pounds et al. 2003; Cappi 2006; Tombesi et al. 2010a b, and 
references therein). In particular. X-ray BALs in He-like iron, 
blueshifted up to Vgut - 0.7c, have been observed in gravita- 
tionally lense d BAL and mini-BAL QSOs (IChartas et al.ll2002l 
[20031120093) . 

Theoretically, powerful winds can be launched from AGN 
accretion disk either by thermal, magnetic, or radiative pressure. 
In either case, a wind will be launched only if the pressure over- 
comes the SMBH gravity pull. As a rule of thumb, the closer to 
the SMBH the launching point is, the higher the wind terminal 
velocity, (see e.g. KonigI ,2Q0.6; Eyerett , 20Q7 ; Proj ga 2007, for a 
review). Observations of the wind outflow velocity are then cru- 
cial in order to constrain the launching radius, and hence the 
physical driving mechanism. The thermal mechanism can ac- 
count for low-velocity outflows arising at large radii (e.g. the 
X-ray warm absorber. Chelouche & Netzeil2005l) . Instead, in or- 
der to explain the extreme velocities observed in the other ab- 
sorbers, either magnetic or radiation pressure must be invoked. 
For instance, line driven accretion disk wind models show that 
at high accretion rates, the UV radiation pressure from the disk 
can be high enough to launch powerful winds from distances of 
a few tens of rg. In the inner part of this flow, high column den- 
sities of X-ray absorbing gas shield the wind from the strong X- 
ray continuu m source, p reventin g it from b ecoming overionized 
dMurrav et al. 1995; Pr oga et al.l i2000; Proga.&Kallman 2004). 
Thus, the geometry and the dynamics of such winds depend crit- 



ically on the UV/X-ray flux ratio and on the AGN accretion rate. 
On the other hand, for the launch and acceleration of magnet- 
ically driven disk winds, there is no need for X-ray shielding, 
and the accretion rate is not a critical parameter; however, the 
ionization state of the gas, which is a critical parameter, wil l still 
depend on them (e.g.. [Konigl & Kartielll994l;lEverettll2005h . 

BAL and mini-BAL QSOs are actually known to be "X-ray 
weak" with respect to the average QSOs. A comparison of the X- 
ray fluxes of BAL and non-BAL QSOs with the same UV flux 
indicates that BAL QSOs typically have X-ray fluxes that are 
10 to 30 times lower th an those of non-BAL QSOs (Green e t al 
'1995; Laoret al. 1997). A numb er of studies (e .g. lBrandt et al 



,2000; Green et al, 2001; Gallagher etaP l2006h strongly sup- 
port absorption as the reason for the X-ray weakness. Recently 
a number of "X-ray bright" BAL and mini-BAL QSOs have 
been discovered, with the latter class X-ray brighter and less X- 
ray absorbed than the former (iGiustini et a l. 2008; Gibson et all 
2009a). The NAL QSOs have been shown to be generally unab- 
sorbed in X- rays, and "X-ray brig hter" than both mini-BAL and 
BAL QSOs (' Chartas et al.ll2009al) . 

The physical link between the UV and X-ray outflowing ab- 
sorbers and the dynamical behavior of the wind are far from be- 
ing understood. Reported estimates of the distances of the ab- 
sorbers from the centra l SMBH range from the inner regions of 
the accretion disk (e.g., Chartas et al. 2009c) to the parsec-scale 
torus (e.g., iBlust in et al. 2005) to the kiloparsec scales (e.g., 
iDunn et a Currently the uncertainties of the launching 

radii are relatively large, which translates to a high uncertainty 
on the mass outflow rate and on the kinetic energy injection 
associated with such winds. It could also be that the launch- 
ing radii depend on the mass of the SMBH. Constraining the 
physical mechanism responsible for launching and accelerating 
AGN winds, i.e. understanding the link between the accretion 
and ejection processes in AGN, can help not only in understand- 
ing the evolution of these sources, but also in quantifying their 
impact on the s urrounding environment, i.e. the amount of feed- 
back (see e .g. iDi Matteo etaP l2005t IScannapieco et alJ 120051; 
[Eivisl[2006h . 

In this article we present the results of an HMM-Newton 
observational campaign on the mini-BAL QSO PG 1126-041, 
which provides both the largest dataset (four pointings) and the 
deepest X-ray exposure (130 ks) on a mini-BAL QSO to date. 
The source is extremey interesting for showing all three ab- 
sorbers discussed above: mini-BALs in the UV, an X-ray warm 
absorber, and a highly ionized, high-velocity X-ray outflow. 
Constraining the X-ray/UV properties of this source can give 
precious insight into the physics of AGN winds. 

The structure of the article is as follows. Section [2] summa- 
rizes the source properties known from the literature. Section [3] 
reports on the new XMM-Newton observations and data re- 
duction. EPIC pn timing analysis results are first presented in 
Section |4| EPIC pn and MOS spectral analysis is then pre- 
sented in Sections [5] (average spectra of each epoch) and [6| 
(time-resolved 2009 Long Look pn spectra). Section [7] contains 
the results of simultaneous optical/X-ray photometric analysis. 
Results are discussed in Section [8] and conclusions presented in 
Section|9| 

A cosmology with Ho = 70 k m s-1 Mpc"', Qa - 0.73 and 
Qm = 0.27 (iKomatsu et al.ll20Tl1) is used throughout the paper 

2. PG 1126-041 

PG 1126-041 (a.k.a. Mrk 1298) is a nearby (z = 0.06), radio- 
quiet AGN with a luminosity in between those that are typ- 
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Table 1. Main properties of PG 1 126-041 



RA 


llh 29m 16.6s 


(a) 


DEC 


-04h 24m 08s 


(a) 


z 


0.062 


(a) 


Dl 


284 Mpc 




K 


4.35 X 10 cm 


(b) 


E(B - V) 


0.055 mag 


(c) 


Mbh 


5.0 - 12 X 10' Mo 


(d-e) 


mg 


14.92 mag 


(a) 


Lhoi 


8.9x 1O"L0 


(f) 


L2-10 


2 X lO''^ erg s"' 


(h) 


IIV 


5000 km S-' 


(g) 




16500 km S-' 


(h) 



References: (a ) ISchrmdt & GreenHl983h: (b)lKalberla et alj ( 120051) : (c) 

Schlegel et al.7 l998^; (d)IVestereaard & PetersonI |l2006l); (e) 

basyra et al. (2007); (f) Surace et al. ( 2001); (g) Wang et 3] ( 119991) ; 
(h) This work, referred to the 2009 Long Look observation. 



ical of Seyf erts and QSOs (Mr = -22.8, Lhoi ~ 1,0'^^» ' 
ISchmidt & Gr een 1983 dSanders et aL 1989; S urace et al.ll2001l) . 



and weak [O III] line 
characteristics of 



emission 
narrow line 



Its optical strong Fe I I 
(iBoroson & GreenI Il992h are 
Seyfert 1 galaxies (NLSl), although the FWHM of its H/3 line 
of 2150 km s"' is sl ightly more than the 2000 km s"' value used 
in defining NLSls dShuder & Osterbrockl[T98l . Its black hole 
mass is estimated to be logiMsn/Mo) ~ 7.7 using the L,i(5 lOOA) 
continuum luminosity density and the FWHM(Hy6) measured 
from a single-ep och optical spectrum an d the empirical rela- 
tions reported by IVestergaard & PetersonI (i200 6). However, the 
stellar velocity dispersion in the host galaxy of PG 1126-041 
measured from the stellar CO absorption line width in the 
VLT H-band spectra gives a slightly hig her black hole mass of 
Iog(MBH/Mo) -8.1 (iDasvra et aLli2007l) . 

PG 1126 - 041 is also classified as "soft X-ray weak'O] by 
iBrandt et al.l (l2000h . who measured an ~ -2.19 based 
on ROSAT X-ray observations and on the optical flux density 
as given by Neugebauer et al. (1987). Tanan baum et al. (198 6) 
measured an qtox - -190 referring to an earlier Einstein X- 
ray observation and the 2500 A rest frame flux density as given 
bylSchmidt & Green ( 1983). Neither of these ffox measurements 
were based on sim ultaneous X-ray and optical/UV observations. 
IWang et alj d 19991) report on the lUE and ROSAT observations of 
PG 1 126-041. The 1992 lUE observation shows a high-velocity 
UV absorbing outflow as mini-BALs in the C IV, N V, and Si IV 
ions blueshifted up to ~ 5 000 km s ' . The simultaneous ROSAT 
pointing revealed a 0.2-2 keV spectrum best-fitted by a steep 
powerlaw emission with F ~ 2.8, absorbed by a column density 
iVw ~ 3 X 10^^ cm"^ of ionized gas, log^ ~ 1.7 erg cm s"'. 
A subsequent lUE spectrum taken in 1995 showed a continuum 
that was brighter by a factor of two and weaker absorption lines. 
iKomossa & Meerschweincheni ([2000) confirm the detection of 
ionized absorption in the 1992 ROSAT spectrum and also sug- 
gest X-ray temporal variability on timescales as short as 800 s. 
The main properties of PG 1 126-041 are summarized in Table[T] 



3. Observations and initial data reduction 

XMM-Newton observed PG 1126-041 four times, once in 2004 
December, twice in 2008 June and December, and once more 



' Soft X-ray weak sources are those being characterized by an ob- 
served optical-to-X-ray spectral index q-qx < -2, see SectionlT] 



in 2009 June. In Table |2] we list the main parameters of each 
dataset. We reduced all the datasets using the XMM-Newton 
SAS V. 10.0.0 and calibration files generated in August 2010. For 
each dataset we extracted the light curve from the whole field 
of view of events with energies greater than 10 keV for each 
European Photon Imaging Camera (EPIC) MOS cameras, and 
with energies 10 keV < £ < 12 keV for the EPIC-pn instrument. 
We then filtered out the time periods in which there was signif- 
icant background flaring. Source counts from the cleaned event 
files were extracted from circular regions centered on the PG 
1 126-041 position, with ~30^5" radii, depending on the signal- 
to-noise ratio (S/N) as determined with the eregionanalyse 
task. Background counts were extracted from circular regions 
of the same area as the source ones. We retained up to double 
events (flagged as #XMMEA_EP) for the pn dataset and quadru- 
ple events (flagged as #XMMEAJ:M) for the MOS datasets. 
For each dataset we generated the ancillary response file and the 
redistribution matrix file at the source position using the tasks 
arfgen and rmfgen. Source light curves were created in dif- 
ferent energy bands and corrected with the epiclccorr task 
for both the local background and the mirrors/detector ineffi- 
ciencies. Due to the low X-ray fluxes involved, pile-up effects 
are negligible in our spectral analysis. For the same reason we 
cannot make use of the reflection grating spectrometer (RGS) 
aboard XMM-Newton. Optical monitor (OM) data were all taken 
in Image mode and were processed using the standard procedure 
as suggested by the XMM-Newton SOC. 

4. Timing analysis 

The lefthand panel of Fig. [1] shows the background-subtracted 
0.2-10 keV pn light curves of PG 1126-041 for each epoch of 
observation, binned to 1 ks. Compared to the light travel time 
f£ = rg/c, the time bin corresponds to ~ 1.5 - 4r^ for the two 
different estimates for the black hole mass of PG 1 126-041 . The 
variations on month timescales are dramatic: the source flux in- 
creases by a factor of ~ 4 from December 2004 to December 
2008, then comes back to a low count-rate regime in June 2009. 
For comparison, the average 0.2-10 keV pn spectra of each 
epoch of observation are shown in the righthand panel of Fig. [T] 
Most of the spectral variability on month time scales occurs at 
energies E < 6 keV. The spectral shape is peculiar and obviously 
deviates from a simple power-law. A prominent broad absorption 
feature is evident in all the spectra around E ~0.6-1.5 keV. 

The high S/N ratio and the long contiguous good exposure 
time allowed us to perform a more detailed timing analysis on 
the 2009 Long Look pn exposure. Figure |2] summarizes the re- 
sults. The lefthand panel shows the root m ean square variabil- 
ity (RMS; Yaughan et al. 2003; .Ponti et al., 2004) computed on 
the 0.2-10 keV band: on very short time scales (ks), the source 
flux varies up to 50% at energies £ > 1 .5 keV, while it is con- 
stant at lower energies. The middle panel shows the background- 
subtracted lightcurves extracted in the 0.2-1.5 and 1.5-10 keV 
bands; following the 1.5-10 keV band flux fluctuations, we split 
the exposure in eight time intervals on which we perform tem- 
porally resolved spectral analysis. Three representative spectra 
extracted during the 2nd, the 7th, and the 8th time intervals are 
shown in the righthand panel, and most of the spectral variability 
on ks time scales occurs at 1 .5 keV < E <9 keV. 

From the timing analysis, it follows that the X-ray spectra of 
PG 1 126-041 can be divided into two spectral components: one 
dominates at £ > 1.5 keV and contributes 100% to the ks time 
scale variability, while the other one dominates at lower energies 
and is significantly variable only on longer (month) time scales. 
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Table 2. X-ray observation log 



OB SID 


Date 


Duration 


Exposure 


Mode 


Filter 


Count-rate 


log/b.2-10 


Epoch 








Ml/M2/pn 






MOS/pn 


pn 








[ks] 


[ks] 






[0.1 ct s-'] 


[erg cm"^ ' ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


0202060201 


12/31/2004 


33.8 


32.9/32.9/28.3 


LW 


t 


0.61±0.02/2.01±0.03 


-11.88 


Dec. 2004 


0556230701 


06/15/2008 


31.4 


3.7/3.5/3.2 


FF 


m 


0.91±0.05/ 3.01 ±0.10 


-11.70 


Jun. 2008 


0556231201 


12/13/2008 


11.9 


4.9/4.7/3.7 


FF 


m 


1.68±0.06/5.93±0.20 


-11.52 


Dec. 2008 


0606150101 


06/21/2009 


134.3 


85.4/89.2/91.2 


FF 


m 


0.57±0.01 / 1.96±0.02 


-11.93 


2009 Long Look 



Notes. Col.(l): Observation ID; Col.(2): Date of observation; Col. (3): Nominal duration of tlie observation; Col. (4): Net exposure time for eacli 
instrument after tlie background flaring filtering was applied; Col. (5): Observing Mode, LW= Large Window, FF=Full Frame; Col. (6): Optical 
Filter applied, t=thin, m=medium; Col. (7): Net count-rate in the 0.2-10 keV and 0.3-8 keV range for the EPIC-pn and EPIC-MOS instruments, 
respectively, after the local background subtraction; Col. (8): Observed EPIC-pn flux in the 0.2-10 keV band; Col. (9): Epoch name as used in the 
text. 
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* 
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Time (s) Time (s) Time (s) 



5x10* 



■ DEC 2004 
• lUN 2008 
▲ DF.C2IIIIN 
JUN 2009 



Time (s) 



0.5 1 2 5 

Observed Frame Energy (keV) 



Fig. 1. Variability of PG 1 126-041 on month time scales. Left panel: background subtracted 0.2-10 keV pn light curves extracted in 
the four different epochs of XMM-Newton observations and binned to 1 ks: December 2004 (blue squares), June 2008 (red circles), 
December 2008 (green triangles), and 2009 Long Look (black circles). Right panel: the time-averaged pn spectra corresponding to 
the time intervals shown on the left and using the same color codes. 




0.2 0.5 1 2 5 10 

Observed Frame Energy (keV) 



■ti 



V 

4 ! 6 7 



5x10* 
Time (s) 




Observed Frame Energy (keV) 



Fig. 2. Variability on kilosecond time scales detected during the 2009 Long Look observation. Left panel: RMS variability computed 
for the pn observation over the 0.2-10 keV band; arrows mark upper limits, error bars are at Icr confidence level. Middle panel: 
0.2-1.5 and 1.5-10 keV background-subtracted pn light curves binned to 2 ks; the numbers in the 1.5-10 keV light curve mark the 
eight time intervals on which time-resolved spectral analysis was performed. Right Panel: three representative spectra extracted 
during the 2nd, the 7th, and the 8th time intervals of the 2009 Long Look observation. 



With these considerations in mind, we proceed with the analy- 
sis of both the average (Sect. |5li and the time-resolved (Sect. |6]l 
spectra. 



5. Spectral analysis: the average spectra 

In our spectral analysis we used the Xspec v. 12.6.0 soft- 
ware. All the spectra were grouped to a minimum number of 
20 counts per energy bin in order to apply the^^ statistics in the 
search for the best-fit model (Avni 1976). Errors are quoted at 
Icr confidence, unless otherwise stated . All the models include 
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3 
O 

U 
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X 



in 
o 
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0.5 1 2 5 
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in 
I 



: \ 1 1 — 1 — r 


1 .^•••^••W^ 


i 1 — 1 — 1—1-4 


-H-l 1 1 1 h- 


1 1 1 r 






\ 1 1 — h- (-4 














fir"' N 



0.5 1 2 5 
Observed Frame Energy (keV) 

Fig. 3. Top panel: average 2009 Long Look pn (black), MOS 1 
(red), and MOS2 (blue) spectra modeled with a power-law emis- 
sion absorbed by the Galactic interstellar medium, along with 
the residuals to the model, scaled by the error bar. Bottom panel: 
spectral residuals to the power-law model for the other three 
epochs: Dec. 2008 (green triangles), Jun. 2008 (red circles), and 
Dec. 2004 (blue squares); the data have been rebinned for visual 
purpose only to respectively 5, 5, and 7cr. 



Galactic absorption b y a column density A^^°' - 4.35x10^" cm - 
dKalberla et al.l l2005l) . We focused our analysis on the 0.2-10 
keV band of the pn data and on the 0.3-8 keV band of the 
MOS data. We simultaneously fit the average spectra of the four 
epochs, Dec. 2004, Jun. 2008, Dec. 2008, and 2009 Long Look, 
to investigate the causes of the observed variability on month- 
long time scales. For each epoch, the model parameters were 
kept fixed between the MOS and the pn datasets, except for 
the primary power-law normalizations that were left free to vary 
within 10% to account for the instrumental cross-calibration un- 
certainties. For the sake of clarity, in most of the following plot 
we only show the pn data. 

The spectra were initially fit to a simple power-law model, 
with fixed slope and free amplitude between different epochs. 
Results are shown in Fig. [3] The fit is very poor, with a chi square 
per degrees of freedom x^/v = 8801/1591, and a very flat pho- 
ton index F ~ 0.6. The spectral residuals of all the four epochs 
have a similar shape, and indeed allowing for the spectral slope 



of the power-law to vary between epochs does not improve the 
fit statistics. Three main deviations from the simple power-law 
model are consistently observed in spectra from all three instru- 
ments and in all four epochs of observation: (a) both negative 
and positive residuals in the soft band, £ < 1 .5 ke V, (b) a deficit 
of counts at the highest energies probed by the EPIC camera, 
£ ~ 7 - 10 keV, and (c) a cutoff' at £ < 2 keV. 

We then tested several models of increasing complexity, to 
reproduce the complex spectral shape of PG 1126-041. A neu- 
tral absorber fully covering the emission source is not required 
by the data. On the other hand, allowing for the absorber to only 
partially cover the source (zpcfabs model in Xspec) signifi- 
cantly improves the fit statistics with respect to the power-law 
model (Ax^/Av = 5148/8), but still gives an unacceptable fit 
Otr ~ 2.3). A better fit for the same number of degrees of free- 
dom (A;^'^/Av = 6567/8) is given by a ionized absorber totally 
covering the source (modeled with XstaiQ), but again, the fit is 
statistically unacceptable (x^ ~ 1.5). 

The ionized absorber was then allowed to only partially 
cover the X-ray source (model A), by adding a secondary soft 
power-law with the same slope as the primary one. The ratio be- 
tween the normalizations of the two power-laws gives the frac- 
tion of the direct (unabsorbed by intrinsic absorption) flux F; 
the absorber covering factor is then C/ = (1 - F). We obtain 
an improvement in the fit statistics of A;^'^/Av = 341/4, and the 
resulting fit is marginally acceptable, with ?l ~ 1-199. The 
model A parameters are reported in Table [3] along with their Icr 
statistical eiTors. We note how the direct flux fraction in every 
epoch is around 2-3%, as is observed for the scattered compo- 
nent in absorbed AGN, and scattered emission accompanied by 
recombination emission lines would be naturally associated to 
the presence of ionized gas along the line of sight. Unfortunately, 
with the S/N of our observations we cannot distinguish among 
the partial covering and scattering scenarios in reproducing the 
soft band residuals, and the secondary soft power-law physical 
interpretation is not unambiguous. To the first order, model A is 
able to reproduce the complex spectral shape of PG 1126-041, 
in particular by recovering the cutoff at £ < 2 keV and most of 
the soft band spectral residuals. While the Jun. 2008 and Dec. 
2008 spectra are generally reproduced well, the low flux states 
of the 2009 Long Look and Dec. 2004 still show strong residuals 
at £■ < 1 keV and £ > 6 keV. 

In particular, the top panel of Fig. |4] shows a magnified 
view of the 4-10 keV residuals of the 2009 Long Look pn data, 
relative to the simple power-law emission (top panel) and to 
model A (lower panel). We make an initial, phenomenological 
model for these residuals by adding Gaussian profiles with neg- 
ative normalization to model A of the 2009 Long Look dataset. 
Even though the model adopted is only a rough approxima- 
tion of a more complex physical context, the basic properties 
such as the FWHM and the EW of the absorption features are 
still physically meaningful. Including two narrow (unresolved, 
cr = Q eV) Gaussian absorption lines improves the fit statis- 
tics by A;^'^/Av = 60/4 with respect to model A. The en- 
ergies and equivalent widths of the two absorption lines are 
£f = 7.02+°°? and Ef' = 7.36!°°2 kgV, EWf' = -80!^° 
eV, in the source rest frame. Allowing the 



and EWf' 



-0.02 



127+24 



We applied the precompiled grid25, publicly available at 
ftp : //legacy . gsfc .nasa . gov/sof tware/plasma_codes/xstar/xspect 
which is computed for an ionizing continuum with F = 2, a gas shell 
with n = 10'" cm"^, a turbulent velocity v , urb = 200 km s"', and 
solar abundances, see iKallman & Bautistal ( 120011) and the Xstar 
documentation. 
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Table 3. Spectral analysis results of the simultaneous fit to the 
average spectra of each epoch 



Rest Frame Energy (keV) 
4.24 5.3 6.36 7.42 8.48 9.54 10.6 





Jun. 2009 


Dec. 2008 


Jun. 2008 


Dpr 9004 




Model A: ionized partial covering 






Xspec : 


zpo2 + XSTAR*zpoi 




r 




2.15 


+0.01 
-0.01 




^IkcV 




1 89+''0'' 


1 40+""^ 

^■^"-0.06 


1 00+0°^ 




14 36+° *' 


4 01+0.24 


7 Q7+0.36 
'■^'-0.24 


9 78+0 '** 

^■'"-0.08 


log^"' 


1 66+°"^ 


1 44+0.06 


1 co+0.02 

l.J0_QQ2 


1 C7+O.OI 
^■-"-0.01 




97+0 '" 

"•"-0.01 


98+0 01 


98+°"' 


97+o"i 
"•^'-0.01 






1893/1579 




Model B: ionized partial covering + highly ionized absorber 




Xspec: zpo2 


+ XSTARl*XSTAR2*zpOi 


r 




2.14 


+0.01 
-0.03 




^IkcV 


1 1 9+0.06 
^•^^-0.13 




1 41+0.09 
^■^^-0.08 


1 01+0 03 


"w 


14.82t« 


T 1 7+0.24 
J.i /„ j2 


6 93+° " 
"■'^-'-i.oo 


q 7 c +0.24 
^•'-^-0.10 


log^""'- 


1.66!°:«| 


1 9q+0.04 
'-■^'^-0.04 


1 ^4+0.02 
^■-^^-0.04 


1 C7+0.01 

^•-^ -0.01 






98+""' 
"■^°-o.oi 


98+°"' 
"■'^"-o.oi 


97+°"' 

"•=^'-0.01 




75. 18^^^ 


< 6.82 


< 6.15 


62 54+" ''' 
"^••^ -35.44 


log^'" 


T 4^+0.09 
-'•^-'-0.08 


< 3.00 


< 2.73 


T 7C+O.I2 

'-^-0.20 




054+0"' 


O.O'' 


O.QP 


0.(f 






1771/1570 





Notes. Units are the following: [10"^ ph s"' keV"' cm"-] for the 
power-law normalizations A'lkcv; [10^^ cm"^] for the column densities 
Nw,N'^''', and N'^'", [erg cm s"'] for the ionization parameters 
and ; c for the highly ionized absorber outflow velocity Vg,„; the 
photon index T and the covering fraction Cf are adimensional parame- 
ters. Errors quoted are at Icr confidence level, a superscripted F denotes 
a fixed parameter. 



widths of the absorption lines to be free parameters does not im- 
prove the fit (Ax^/Av = 2/2). In the bottom panel of Fig.|4]we 
show the 99% confidence contours for the absorption lines cen- 
troid energy and (negative) normalization, along with the ener- 
gies of the strongest transitions of Fe XXV and Fe XXVI. There 
is a remarkable correspondence between the detected centroid 
energies and the Fe XXV Hea and Fe XXVI Lya transitions, 
both blueshifted by v„,„ = 0.055c, so we identify the two ab- 
sorption lines with these two transitions. 

The moderately ionized absorber (m.i.) responsible for the 
opacity in the soft band has an ionization state too low to ac- 
count for the strong absorption in the iron K band. So we added 
another layer of highly ionized (h.i.) gas modeled again with 
XSTAR. We used the same input parameters as the grid 25, but 
we find a turbulence velocity v^m-jj — 1500 km s 'to better re- 
produce the residuals. Only in the case of the 2009 Long Look 
data we allowed the velocity shift of the highly ionized absorber 
to be a free parameter; given the much lower S/N of the three 
other epochs, we fixed the velocity shift of this component to 
zero. The addition of the highly ionized absorber significantly 
improves the fit with respect to the moderately ionized partially 
covering absorber model (A^^/Av = 122/9), and gives a good 
overall representation of the data (xj. ~ 1.128). Spectral param- 
eters along with their errors are reported in Table [3] model B. 
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Fig. 4. Top panel: 4-10 keV spectral residuals of the 2009 Long 
Look pn data, relative to the simple power-law model (top panel) 
and to model A (bottom panel). The two arrows mark strong 
negative residuals that are present in both cases. Data have been 
rebinned to 3cr significance. Bottom panel: 99% confidence con- 
tours for the rest frame normalizations and centroid energy of 
the absorption lines used to model the residuals; the black con- 
tours refer to the addition of the absorption lines to the simple 
power-law model, while the red ones to the addition to model 
A. The four vertical green lines mark the theoretical energies of 
the strongest transitions of Fe XXV and Fe XXVI, while the two 
dotted blue lines mark the position of the Fe XXV Hea and Fe 
XXVI Lyo' transitions, both blueshifted by Vo,„ — 0.055c. 



The corresponding theoretical model is plotted in top panel of 
Fig. |5j while the spectral residuals to the model are plotted in 
the bottom panel of Fig. |5]for each epoch. 

Significant residuals at £ < 1 keV relative to model B are 
present only in the 2009 Long Look and Dec. 2004 spectra. 
The residuals could be either due to an improper modelization 
of the absorbers or to an unmodeled component that emerges 
in the soft band when the source is in a low flux state. In any 
case, we know from the timing analysis performed in Sect. |4] 
that in this energy range the source shows no spectral variabil- 
ity on short time scales. We find an improvement in statistics 
of Ax^/Av = 42/2 and Ax^/Av = 23/2 for the addition of two 
naiTow Gaussian emission lines with centroid energies and nor- 
malizations fixed between the epochs. The best-fit rest frame en- 
ergies are = 0.55+° °| keV and £™' = 0.91+° °| keV, while 
the equivalent widths EW™" -10-30 eV and EW'^"" -15-40 
eV, depending on the epoch of observation. The fit statistics af- 
ter including the two soft emission lines is^^/v = 1710/1566. 
The observed energies might correspond to O VII K and Fe L 
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Table 4. 2009 Long Look observation time slices over which 
time-resolved analysis was performed. 
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Fig. 5. Top panel: the solid lines show the power-law emission 
after the passage through two absorbers, one of which is partially 
covering the source (model B). This model was fit to the aver- 
age 2009 Long Look (black), Dec. 2008 (green), Jun. 2008 (red), 
and Dec. 2004 (blue) epoch spectra. The dashed lines represent 
the direct continuum components that are not absorbed. Bottom 
panel: average 2009 Long Look (black points), Dec. 2008 (green 
triangles), Jun. 2008 (red circles), and Dec. 2004 (blue squares) 
spectral residuals relative to model B; the data have been re- 
binned for visual purpose only to respectively 15, 5, 5, and lOo". 



emission features, but the exact identification is prohibitive with 
the pn spectral resolution. Physically, possible contributors to 
the soft emission can be a scattered component accompanied by 
recombination lines from a photoionized gas, thermal emission 
from an underlying starburs t component, o r a mixture of the two 
(e.g. IRAS 13197-1627, Miniutti et al. 2007). Unfortunately, PG 
1 126-041 is not detected by the RGS, and the pn spectral reso- 
lution does not allow the different scenarios to be distinguished, 
thus we do not try to model the soft residuals any further 

As for the high-energy residuals, an excess of counts at ob- 
served energies E < 6 keV is present in the 2009 Long Look 
data, suggesting the presence of complexities in the Fe K band. 
We note that at these energies PG 1126-041 is highly variable 
on kilosecond time scales (Sect. |4]l. The addition of an Fe I Ka 
emission line, which is forced to be narrow and has a centroid 
energy of 6.4 keV in the source rest frame improves the fit statis- 
tics by Ax^/Av = 15/1, corresponding to a significant detec- 



Interval 


Exposure 


Count-rate 


/2-10 


/0.2-IO 


1 


21.1 


0.180+0.003 




1 00+""'^ 

-Q Q3 


2 
3 
4 


30.0 
5.4 
7.5 


0.156+0.003 
0.212+0.007 
0.195+0.005 


0^8:8? 

"■"--0.03 
1 17+0.04 
^ • ^ ' -0.05 

98+° '"' 
-Q-!!4 


0.83!«| 

1 ^^+0.04 
^•-'^-0.06 

J 12+""'' 
1 32+^-8^ 


5 


4.6 


0.215+0.007 


1 17+0.04 

-0.05 


6 


5.4 


0.194+0.006 


98+" '''' 


^ ]^ 2+^-8^ 


7 


7.7 


0.338+0.007 




2 31+8:85 


8 


9.2 


0.245+0.006 


1 aQ+H;?^ 


154+8:?^ 



Notes. Col.(l): Time interval as marked in Fig. |2j Col. (2): Exposure 
time in ks; Col.(3): 0.2-10 keV pn count rate in ct s"'; Col.(4): best-fit 
pn observed flux computed over the 2-10 keV observed band, in units 
of 10"'^ erg cm"^ s"'; Col.(5): best-fit pn observed flux computed over 
the 0.2-10 keV observed band, in units of 10"'' erg cm"^ s"'. Errors 



quoted at Icr confidence level. 



tion (F-test probability > 99.97%). The best-fit model statistics 
sae x^/^ - 1691/1565, with a reduced chi square = 1.080. 
There are still shallow residuals in the Fe K band, which could 
be modeled with a broad Gaussian emission line, a relativistic 
Fe K line, or a partial covering absorber, in all cases suggesting 
the possible presence of further complexities in the iron K band 
(see M. Giustini Ph.D. thesis for details). However, given the low 
statistical significance of the residuals, their appearance only in 
the 2009 Long Look observation, and the high variability of the 
source at those energies on a short time scale during the same 
observation, we do not try to model these residuals and turn our 
attention on the time-resolved spectral analysis of the 2009 Long 
Look data. 



6. Time-resolved spectral analysis 

We performed time-resolved spectral analysis on the EPIC pn 
data extracted in the eight time intervals of the 2009 Long Look 
observation marked in the middle panel of Fig. |2] Exposure 
times and background-corrected count rates are reported in 
Table|4]for each interval. Given that the count rates are the same 
in several intervals, we merged the third and the fifth intervals 
together and the fourth and the sixth one together The unfolded 
0.2-10 keV pn spectrum is plotted in Fig. |6]for each of the six 
resulting intervals. 

At energies £ > 1 .5 ke V there are variations of about 50% on 
a few kiloseconds timescales (e.g., from the second to the third 
time interval, and from the seventh to the eighth time interval), 
and variations as high as 100% with respect to the average are 
observed during the seventh interval. All the spectra were simul- 
taneously fitted to the best-fit model B described in the previous 
section, with all the parameters tied between the time intervals 
except for the power-law normalizations. The fit to the data is 
acceptable, x^/'^ - 925/889, however the model fails to prop- 
erly reproduce the Fe K band absorption residuals. Allowing for 
the highly ionized absorber column density and ionization pa- 
rameter to vary within each temporal slice improves the fit by 
Ax^/Av = 45/6. Because of the much lower S/N in individual 
temporal slices, it is not possible to assess whether the absorber 
varied in column density, ionization state, or a combination of 
them. Fluxes computed using this best-fit model are reported in 
Table m 
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Fig. 6. The unfolded 0.2-10 keV EPIC-pn spectra extracted in 
the different time intervals marked in Fig. |2l 1*' (black), 2""' 
(red), 3''' + 5"' (green), 4"' + 6"' (blue), 7"' (light blue), and 8"' 
(magenta). Data points have been visually rebinned to 10 cr con- 
fidence level. 



We then removed the highly ionized absorber from the model 
and visually inspected the residuals in each temporal slice. Given 
the similar shape of the residuals, we merged the third, fourth, 
fifth, and sixth time intervals together, and the seventh and eighth 
time intervals together. The 5-9 keV spectral residuals are shown 
in the left column of Fig. |7]for each of the four resulting inter- 
vals; time is flowing from top to bottom. There is clear variability 
of the highly ionized absorber on very short time scales: in par- 
ticular, there are no signs of high-energy absorption during the 
second interval, while the effect of iron absorption during the last 
time interval is dramatic. The right column of Fig. |7] shows the 
corresponding 68% and 90% confidence contours for the energy 
and (negative) normalization of a narrow Gaussian absorption 
line used to model the residuals. One can see how the complex 
absorption features at ~ 7 keV in the observer frame do disap- 
pear within the few ks that elapsed between the first and the sec- 
ond time intervals, and then re-develop during the third to sixth 
temporal slices. The nearly 7 keV absorption trough becomes the 
deepest and the largest during the seventh-eighth interval, and its 
shape becomes quite complex, as confirmed by the confidence 
contours. The intense absorption structure at ~ 6.5 keV does de- 
velop in a few ks and is also visible in the merged 2009 Long 
Look spectra (top panel of Fig.|4]i. Summarizing, the continuum 
power-law amplitude oscillations dominate the observed spec- 
tral variability on ks time scales. The highly ionized absorber 
is also observed to be variable, while all the other components 
(i.e. the photon index, the partial covering moderately ionized 
absorber, and the naiTow Gaussian emission lines) are found to 
be constant during the observation. 

As for the other epochs, the Jun. and Dec. 2008 observations 
were too short to search for short-term variability. Strong spec- 
tral variability was found instead during the 2004 observation: 
Fig. |8] shows the unfolded pn spectra extracted during the first 
10 ks and during the last 10 ks of the Dec. 2004 observation. 
The variability pattern is the same as of the 2009 Long Look 
observation i.e. the source is variable on ks time scales only at 
energies E > I keV, and its spectral variabilty is dominated by 
continuum amplitude variations. 
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Fig. 7. Left panels: 5-9 keV residuals for each temporal slice 
to model B, after the variable highly ionized outflowing ab- 
sorber was removed. From top to bottom: V (black), 2""' (red), 
6"' (blue), 7"' + 8"' (magenta) interval. Right 



3'-^ + 4"' + 5" 



panels: 68% (continuous line) and 90% (dotted line) confidence 
contours for the energy and intensity of a narrow Gaussian ab- 
sorption line used to model the residuals in each slice. 
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Fig. 8. The unfolded 0.2-10 keV EPIC-pn spectra extracted dur- 
ing the first 10 ks (black) and during the last 10 ks (red) of the 
Dec. 2004 observation. Data points have been visually rebinned 
to 10 cr confidence level. 



7. Optical/X-ray photometry 

Thanks to the presence of the OM aboard XMM-Newton, simul- 
taneous optical/X-ray photometry was performed for each epoch 
of observation. The OM data were all taken in image mode, and 
were reduced and calibrated using standard pipeline procedures 
through the omichain SAS task. Background-coiTected count 
rates were converted in magnitudes for each exposure (exposure 
times are ranging between 1 and 4 ks). We were interested in 
measuring /2500 to compute the optical-to-X- ray spectral index 
ffox = -0.384 log(^25{)0A/^2fcv), as defined bv lTananbaum et al.l 
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(Il979h . We used the weighted mean of the magnitudes mea- 
sured in the filters with effective wavelength straddling 2650 A 
in the observed frame (2500 A in the rest frame). The relevant 
filters are the UVWl and UVM2 (A,// = 2910 and 2310 A, re- 
spectively). We then corrected this magnitude for Galactic ex- 
tinction A^fi-jn = .382, estimated using the extinction law of 
ICardelh et al.l (1198 9). Flux densities for each epoch were then 
computed using Vega as a calibrator Optical photometry results 
are reported in Table |5] 

As for the 2 keV flux density measurements, we used our 
best-fit model to measure the Galactic absorption corrected rest 
frame f2kev for each observation. Flux densities were then con- 
verted to luminosity densities, and the optical-to-X-ray spec- 
tral index was computed. It is well known that the ffox spec- 
tral index in AGN stro ngly depends on the intrinsic UV con- 
tinuum luminosity ( e.g., Vignah et al J2003I: IStrateva et alj2005l 
ISteffen et al.1 120061: lYoung et al.ll2010h . We compared the ob- 
served Oox with the value Q'ox(^2500a) expected on the basis 
of the 2500 A UV luminosity of PG 1126-041, using Eq. 2 
of lStefl^enetaT] (l2006l) . which gives us the parameter A^ox = 
ttox - Q'ox(^2500a)- '■h^'- i'^ '■h^ luminosity range of 

PG 1 126-041 there is a 3cr dispersion in the expected Q'ox(^2500a) 
of about 0. 1, so th at only values of IAq-qxI > 0.1 are statistically 
significant (IStefifen e t al. ,2006). We also computed the f2t'Jv 
minosity densities corrected for intrinsic X-ray absorption, and 
the corresponding aj;™ and Aa'^'^. In Table |6] we report the 
Optical/X-ray photometry results. 

There are variations in the observed luminosity density of 
PG 1126-041 both at UV and at X-ray wavelength. While the 
observed 2500 A flux density steadily increases from 2004 to 
2009, and is varying by about 30%, the observed 2 keV flux 
density increases from 2004 to 2008, and then decreases again 
in 2009, with variations as high as a factor of ~ 4 - 5 between 
the Dec. 2008 and the 2009 Long Look observations. As a re- 
sult, the observed ffox is variable between the different epochs, 
spanning from values typical of radio quiet type 1 AGN (i.e., 
Qox ~ -1.7 in the Dec. 2008 observation) to values typical of 
Soft X-ray Weak AGN (i.e., qtox ~ -2 in the 2009 Long Look 
observation). However, compared to the expected q;ox(^2500a)' 
the source is found to be "X-ray weak" in all the epochs, with 
Affox ~ -0.3 - 0.6. Once the effect of intrinsic X-ray absorp- 
tion is taken into account, the maximum observed variations in 
the 2 keV flux density between the different epochs decrease to 
about a factor of two. Consequently, also the Affox decreases to 
Aa™" ~ 0.1 - 0.3. It is worth noting that there could be intrin- 
sic UV absorption at the redshift of the source, e.g. by dust in 
the AGN host galaxy; the amount of intrinsic UV absorption is 
unknown so we do not try to model it, but we keep in mind that 
this effect might make the actual (intrinsic) aox flatter than Co"". 



8. Discussion 

PG 1 126-041 is observed to be highly variable in X-rays on time 
scales of both months and kiloseconds. The variability on time 
scales of months (observed to be as high as 4x in flux) is domi- 
nated by a spectral component emerging at £ < 6 keV (Fig.[T]). 
On the other hand, the RMS variability analysis performed on 
the 2009 Long Look pn observation (92 ks of contiguous good 
exposure time) revealed the presence of a spectral component 
emerging only at E > 1.5 keV, which dominates the variability 
on time scales of kiloseconds (up to 50% variations in flux. Fig. 

m. 



A wealth of information has been obtained both on the intrin- 
sic X-ray continuum emission of PG 1 126-041 and on the repro- 
cessing media that happen to be in the inner regions of this AGN. 
We discuss our results starting from the inner regions around 
the SMBH (i.e., the shortest time scales) and moving farther out 
(i.e., the longest time scales). We adopt t he SMBH mass esti - 
mate of Mbh = L2 x 10** Mo as given by Dasvra et al. (200?!). 
The Mbh estimated by Vestergaard& Peterson (2006) through 
scaling relations is a factor o f 2.4 lower. Adoptin g the bolomet- 
ric luminosity Ltoi - 10'^ Lq dSanders et alJl989l) and assuming 
an accretion efficiency 77 = 0. 1 gives an accretion rate needed to 
power PG 1126-041 of M ~ 0.7Mo/yr and an Eddington ratio 
A = Lhoi/LEdci ~ 0.26. The gravitational radius is rg ~ 1.8 x lO'-' 
cm and the corresponding light travel time is ti ~ 600 s. 

Two absorbers are detected in the X-ray spectra of PG 1 126- 
041, a moderately ionized (log^'" ' ~ 1.5 erg cm s"') and a 
highly ionized (log^'" ~ 3.5 erg cm s"') one. Strong varia- 
tions in both the intrinsic continuum and of the highly ionized 
outflowing absorber are responsible of the observed kilosecond 
time scale variability. The power-law photon index F ~ 2 is 
found not to vary on kilosecond time scales, while the intensity 
of the intrinsic power-law spectrum ifollows the pattern of the 
1.5-10 keV count rate over the whole observation (see the mid- 
dle bottom panel of Fig. |2] and Table |4]i, and doubles in a time 
interval At lasting about 8 ks (the seventh interval). The variabil- 
ity time scale of the intrinsic continuum can set a constraint on 
the geometrical size D of the X-ray emission region through the 
causality argument, D < cAt ~ 2.4 x 10'"* cm ~ I3rg. The highly 
ionized absorber is observed to be outflowing at Voui ~ 16,500 
km s ' during the average 2009 Long Look observation. In the 
time-resolved spectral analysis of the same exposure, the out- 
flowing absorber is found to be variable over very short time 
scales (Fig.|7]i. With the much lower S/N of each temporal slice, 
it was not possible to assess whether the absober varied in ioniza- 
tion state, column density, or blueshift. However, the very short 
time scale variability suggests that the absorber is very compact 
and very close to the X-ray source and that we are possibly ob- 
serving rapid mass ejections from the inner regions of the accre- 
tion disk, e.g., the base of an accretion disk wind. We note in 
particular how the absorption features do develop in a very short 
time scale, from the sixth to the seventh interval, together with 
the strong continuum flare. 

By making some assumptions, we can try to estimate the 
mass outflow rate Mom associated to the highly ionized absorber. 
In general, the mass outflow rate can be written as 

Moui ocA(r)p(r) v„Jr) 

where A{r) is a geometrical factor that accounts for how the flow 
diverges, and p(r) and Voutif) are the density and velocity of 
the flow. By assuming a spherically symmetric, isotropic, steady 
flow with a constant velocity, one can write 

Mout ^Anr^n nin Vou, C/Fy ( 1 ) 

where niu is the hydrogen atom mass, r the absorber distance 
from the central SMBH, n the absorber number density, C/ = 
dQ.IAn is the solid angle occupied by the flow as seen by the 
central point source, and the volume filling factor 

If we use appropriate values for PG 1 126-041, we can rear- 
range the expression and write 

gj^vMoyr- 



M„„, = 0.33(nii)(rio)^^__, . 

,16,500kms , 

where n\\ - n/lO" cm"^, r\Q - r/lQrg. In the expression above, 
the geometrical covering factor < C/ < 1 and volume filling 
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Table 5. Optical photometry 



Epoch 


V 


U 


UVWl 


UVM2 


UVW2 


M(2500) 




[5430 A] 


[3440 A] 


[2910 A] 


[23 10 A] 


[2120 A] 


[2500 A] 


Dec. 2004 


14.71±0.01 


14.09±0.01 


13.87±0.01 


14.13±0.02 


14.01±0.04 


13.60 


Jun. 2008 


(...) 


(...) 


13.71±0.01 


13.92±0.01 


(...) 


13.42 


Dec. 2008 


(...) 


(...) 


13.82±0.01 


13.86±0.02 


(...) 


13.46 


2009 Long Look 


(...) 


(...) 


13.564±0.004 


13.656±0.008 


13.72±0.02 


13.22 



Notes. For each filter we list the eff'ective wavelength; the magnitudes measured in each filter are not corrected for the Galactic extinction, while 
the extrapolated magnitude at 2500 A is corrected for the Galactic extinction. 



Table 6. Optical/X-ray photometry 



Epoch 


log /2500 


log flkeV 


log ^2500 


log^2feV 






log^2T.V 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Dec. 2004 


-25.46 


-30.46 


29.49 


24.49 


-1.92 


-0.52 


25.28 


-1.68 


-0.28 


Jun. 2008 


-25.39 


-30.17 


29.56 


24.77 


-1.84 


-0.43 


25.31 


-1.63 


-0.22 


Dec. 2008 


-25.40 


-29.80 


29.55 


25.15 


-1.69 


-0.28 


25.57 


-1.53 


-0.12 


2009 Long Look 


-25.30 


-30.46 


29.64 


24.49 


-1.98 


-0.56 


25.26 


-1.62 


-0.20 



Notes. Col.(l): Epoch of observation. Col.(2): 2500 A rest frame flux density corrected for Galactic extinction. Col.(3): 2 keV rest frame flux 
density corrected for Galactic absorption. Col. (4): 2500 A rest frame luminosity density corrected for Galactic extinction. Col. (5): 2 keV rest frame 
luminosity density corrected for Galactic absorption. Col. (6): Observed optical-to-X-ray spectral index. Col. (7): Dilference between the observed 
Q-ox and the one expected on the basis of the ^2500 luminosity. Col. (8): 2 keV rest frame luminosity density corrected for Galactic and intrinsic 
absorption. Col. (9): Optical-to-X-ray spectral index corrected for the intrinsic X-ray absorption. Col. (10): Difference between the corrected or™" 
and the one expected on the basis of the ^2500 luminosity. Flux densities are in [erg s"' cm"^ Hz"'], luminosity densities are in [erg s"' Hz"']. 



factor Fv < 1 are unknown, and depend on the wind geometry 
and duty cycle. The density and the radial distance of the flow 
are also unknown. 

Assuming that the outflowing absorber is a thin spherical 
shell of gas in photoionization equilibrium, one can use the def- 
inition of the ionization parameter ^ = Lion/nr^ to replace the 
unknown density and distance with the ratio Li„„/^. However, in 
addition to the strong assumptions made regarding the spherical 
geometry of the outflow, there are uncertainties in the ionizing 
continuum luminosity, and especially in the actual value of the 
ionization parameter that may affect our mass outflow rate esti- 
mates. 

In our spectral analysis results we quoted only the statistical 
errors on ^: the systematic errors related to the input parame- 
ters used in the photoionization code (Xstar in our case) can be 
much larger, especially for highly ionized absorbers where only 
a few spectral transitions are available (e.g., Fe XXV, Fe XXVI). 
The major dependence of the best-fit ^ value is related to the 
slope r,o„ adopted for the ionizing power-law continuum. This 
is because ^ oc ny/nu, where is the ionizing photons num- 
ber density, rtu the photoionized gas number density. By fixing 
the ionizing luminosity, a steeper r,v,„ results in a smaller num- 
ber of high-energy photons able to strip electrons from highly 
ionized Fe than for a flatter r,o„. This leads to higher best-fit val- 
ues measured for ^ when using a steeper ionizing continuum. 
To show this dependence, we generated several Xstar absorber 
models using the same input for all the physical parameters, but 
different slopes F,,,,, . We then fitted the 2009 Long Look average 
spectra to each of these models, recording the best-fit ^ value 
measured for the highly ionized absorber responsible for the ab- 
sorption features visible in Fig.H] Results are plotted in Fig.|9l 
the systematic error on ^ is very large, and even restricting the 
slope of the ionizing continuum to values typical of type 1 AGN, 



2.2 2.4 



Fig. 9. Measured ^ for the highly ionized outflowing absorber 
in the 2009 Long Look spectra, with different F,o„ used as in- 
put in the Xstar simulations; Error bars are purely statistical 
errors at Icr confidence level. The vertical dotted line marks the 
best-fit measured F for the 2009 Long Look spectra, along with 
its statistical (dashed lines) and systematic (continuous line) er- 
rors that have been estimated considering the uncertainties of the 
power-law normalization and of the ionized absorbers. 



1.8 < F,o„ < 2.2, gives a systematic uncertainty on ^ of a factor 
of about 30. Using the measured slope F for PG 1 126-041 along 
with its systematic error gives an uncertainty on our measured 
^h.i. q£ ^ factor of ~ 6. 

Just as for a comparison with other literature results, we can 
roughly estimate the kinetic luminosity La,„ associated to the 
PG 1126-041 outflow as follows. We take the best estimate for 

from Fig.|9] considering the best-fit measured photon index 
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r = 2.15 as the slope of the ionizing continuum. For the ionizing 
luminosity, we take the measured unabsorbed 2- 10 keV luminos- 
ity listed in Table[T] Substituting the ratio of these two quantities 
to the product r^n in Eq.[l] and using the best-fit outflow veloc- 
ity, one can write 

Uin - \mou,vL, ~ 3 (^j [y^] X 10« erg s"', 

i.e., less than one hundreth of the bolometric luminosity of PG 
1126-041. 

The uncertainties on the physical quantities related to the 
outflow are huge, and the assumptions made to derive the ex- 
pression for the mass outflow rate and kinetic luminosity are 
likely to be incorrect. The presence of the accretion disk breaks 
the spherical symmetry by itself; furthermore, the highly ionized 
absorber is observed to be variable on very short time scales: the 
geometrical and dynamical effects are most probably very im- 
portant in accretion disk wind scenarios. For these reasons, it is 
very difficult to give a meaningful (i.e. with less than two or- 
der of magnitudes uncertainties) estimate of the mass outflow 
rate and of its kinetic efficiency. Although one can not exactly 
assess the impact of the fast outflow of PG 1 126- 041 on its en- 
vironment, we stress that its detection is statistically solid and 
model-independent. 

Turning our attention to the larger scales, we note that at 
£■ < 1.5 keV the power-law variability is smeared out by the 
presence of a soft constant component, and at £ < 1 keV the 
highly ionized absorber opacity drops (see top panel of Fig.|5]), 
so explaining the constant soft flux measured during the whole 
2009 Long Look observation. As for the soft component, we 
modeled it as the direct continuum emission that escapes un- 
absorbed by a partially covering ionized gas. However, the low 
fraction of unabsorbed flux, (1 - C/) ~ 0.03, could also physi- 
cally correspond to a scattered component, that would be natu- 
rally associated to the presence of ionized gas along the line of 
sight. With the moderate pn spectral resolution it is not possible 
to distinguish among the different scenarios, and unfortunately 
PG 1126-041 is not detected by the RGS. However, we stress 
that the adopted modeling for the soft band constant component 
influences neither the absorbers nor the continuum physical pa- 
rameters. 

The PG 1126-041 X-ray variability on long time scales is 
dramatic, and is especially pronounced at energies E < 6 keV. 
The /2kev is observed to be variable up to 4.5 higher in the 
six months elapsed between the Dec. 2008 and the 2009 Long 
Look observations (Fig.[T]i. The moderately ionized absorber is 
detected in every XMM-Newton observation of PG 1126-041. 
Its column density varies significantly between the minimum 
N'^^^ ~ 3 X 10^^ cm"^ measured in the Dec. 2008 observation 
and the maximum N'^'' ~ 1.4 x 10^^ cm"^ measured in the 
2009 Long Look observations. In fact, the moderately ionized 
absorber column density variations dominate the spectral vari- 
ability observed on month time scales, together with amplitude 
variations of the intrinsic continuum up to a factor of two (Table 
[3]). The ionization state of the moderately ionized absorber 
log^"'-'- ~ 1.3 - 1.7 is typical of the X-ray warm absorbers, how- 
ever the measured column density are much high er than what 
is usu ally observed in warm absorbers (see e.g. iMcKernan et aP 
l2007h . 

A possible physical scenario is that of a UV line-dri ven ac- 
cretion disk wind like t he one developed by Proga & KallmanI 
( |2004 and modeled by ISim et all (|2010|) . These authors com- 
puted the simulated X-ray spectra for a nonspherical, hydrody- 



namical accretion disk wind in an AGN with Mbh — 10 Mq, ac- 
creting with A - LILEdd - 0.5, i.e., with conditions that are not 
too different to those observed in PG 1126-041. The simulated 
spectra change significantly for different inclination angles and 
different snapshots of the flow, confirming both the nonspherical 
character and the highly dynamical behavior of the X-ray spec- 
trum associated with such winds. The spectral signature of such 
flows includes both absorption and emission features. For exam- 
ple, the iron K band is shaped both by absorption by a highly 
ionized (Fe XXV, Fe XXVI) outflowing phase of the wind, and 
by emission and scattering off the highly ionized base of the flow, 
which produces a Compton bump that mimics a red-skewed iron 
K line. In the soft band {E < \ keV), there are contributions by 
a blend of emission lines from elements such as C, O, Fe, and 
bremsstrahlung emission by the wind. 

From a dynamical point of view, the same model also pre- 
dicts a flow component that qualitatively shows ffie same be- 
havior of the highly ionized absorber. At intermediate inclina- 
tion angles above the disk plane and above the "skin" of the 
fast wind, where the gas is most exposed to the X-ray ioniz- 
ing source, recurrent instabilities with blobs (or "puffs") of high 
density can develop. Some of these puffs quickly become overi- 
onized and fail to become part of the wind, while others can be 
able to be efficiently accelerated to velocities even higher than 
those of the fast wind and well above the local escape velocity. 
This flow component would be observationally detectable in the 
X-ray spectra via large column density variations and sporadic 
high-velocity ejections, on short timescales, similar to what was 
observed during the 2009 Long Look observation of PG 1126- 
041. This may suggest that radiation-driven accretion disk wind 
models can account in a self-consistent way for most of the ob- 
served X-ray spectral features in the mini-BAL QSO PG 1126- 
041. Models for the X-ray spectra of accretion disk winds that 
can be fitted to real data are still under developement, there- 
fore we must confine ourselves to qualitative considerations. 
However, we note that in this scenario the radial flow assumption 
is invalid, and the mass outflow rate associated with the highly 
ionized phase is overall very different than in the spherical case. 

The moderately ionized absorber, on the other hand, is not 
variable on short time scales but is strongly variable over time 
scales of months. Its observed variability is consistent with the 
scenario suggested to explain the variability in the UV absorp- 
tion lines of the mini-BAL QSO HS 1603-H3820 over rest fra me 
time scales of weeks and months ( iMisawa et al .l2007bll20Toh . In 
this scenario the variability of the UV absorption lines is the re- 
sult of fluctuations in the continuum, which is caused by a screen 
of clumpy, highly-ionized gas between the UV absorber and the 
continuum source. This clumpy screen has similar properties as 
the moderately ionized absorber we have found in PG 1 126-041 . 

The observed optical-to-X-ray spectral index is found to be 
highly variable, following the moderately ionized absorber vari- 
ability. Although Oox never gets < -2, when compared with the 
expected a'o.v(^2500A) b^sed on the UV continuum luminosity, it 
is found that PG 1 126-041 is observed to be "soft X-ray weak", 
Aaox = -(0.3 - 0.6) depending on the epoch of observation. 
After correcting for the intrinsic X-ray absorption, the source is 
still sUghtly X-ray weak, ^aZ"' = -(0.1 - 0.3). The observed 
Qox variabi lity driven by X-r ay absorption variability is simi- 
lar to what iBallo et al.l (|2008|) observed in the mini-BAL QSO 
PG 1535-1-547. The 2-10 keV unabsorbed luminosity is still very 
low, (L2-10) ~ 2 X 10'*"' erg s"', and compared with Ltoi gives a 
very high bolometric correction, ki,oi ~ 200, that is however sti ll 
compatible with the Kboi - a ox relation as in lLusso et al ] d2010h . 
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9. Conclusions 

Our XMM-Newton observational campaign on the mini-BAL 
QSO PG 1 126-041 allowed us to characterize the complex high- 
energy spectral behavior of the source with unprecedented sen- 
sitivity. We analyzed data for four different pointed observations 
that in all span 4.5 yr. 

The most evident and model-independent result is the detec- 
tion of high column densities of ionized gas along the hne of 
sight, which contribute to the observed spectral variabihty of the 
highly ionized and moderately ionized phase on kilosecond and 
month time scales, respectively. In particular, the highly ionized 
absorber is found to be variable on very short (a few ks) time 
scales, suggesting an origin very close to the SMBH. Also the 
intrinsic X-ray continuum emission is observed to be variable in 
intensity, on both short and long time scales. 

Overall our findings are quahtatively consistent with 
radiation-driven, accretion disk wind model predictions, where 
one expects a high column density of X-ray absorbing gas 
shielding the portion of the wind that is accelerated by UV pho- 
tons. The hot X-ray component of the wind is also expected to 
be highly variable with time in such a model. The mini-BAL 
appearance of the wind of PG 1126-041 is quahtatively consis- 
tent with an orientation eff'ect, in which our line of sight and the 
plane of the accretion disk make a larger angle than when view- 
ing classical BAL QSOs. We are looking at a narrower range of 
UV velocities through the wind (hence the small width of the 
UV absorption features) and at the transition zone between the 
fast UV wind and the hot X-ray component of the flow, where 
we expect strong X-ray spectral variability. 

The present X-ray observational campaign of the mini-BAL 
QSO PG 1126-041 has demonstrated that long and short term 
variability studies are quite powerful tools for identifying the 
physical mechanisms at work and for mapping the dynamics of 
the inner accretion/ejection flow in AGN. 
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